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SUMMARY 


A special  blsignal  (pressure-temperature)  transducer  has  been 
developed  which  measures  the  contact  zone  temperature  between  a pair  of 
rolling  steel  disks.  Temperature  data  have  been  recorded  for  the  follow- 
ing four  lubricants: 

(1)  Synthetic  paraffinic  lubricant 

(2)  Traction  fluid 

(3)  Polyphenyl  ether 

(4)  MIL-L-23699  lubricant. 

The  conditions  for  the  experiments  include  two  loads  (.7  and  1.1  GPa) , 
two  speeds,  and  several  slip  conditions. 

Very  large  temperature  rises  have  been  seen  in  many  of  the  experi- 
ments. For  example,  under  pure  rolling  conditions,  temperature  in  excess 
of  30  C rise  have  been  measured.  For  high  sliding  conditions,  the  tempera- 
ture may  rise  to  200  C. 
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INTRODUCTION 


The  purpose  of  the  research  project  has  been  to  measure,  evaluate, 
and  analyze  temperatures  under  elastohydrodynamic  (EHD)  contact  conditions. 

EHD  contact  occurs  in  many  types  of  rolling  or  sliding  contact  situations, 
such  as  roller  or  ball-bearing  contacts,  gear  tooth  meshes,  and  rolling  disks. 
The  temperature  measurement  technique  developed  during  the  project  involves 
the  use  of  a titanium  temperature  transducer  (thermistor)  coated  on  one  of 
a pair  of  disks  rolling  in  lubricated  contact.  The  change  in  resistance 
of  the  titanium,  due  to  temperature,  is  detected  for  -ach  passage  of  the 
transducer  through  the  contact  region. 

In  addition  to  the  temperature  measurements,  measurements  of  pres- 
sure were  made  using  a manganin  pressure  transducer.  The  temperature  and 
pressure  transducer  were  located  on  the  disk  in  such  a manner  that  the  con- 
tact pressures  and  temperatures  could  Be  detected  simultaneously.  With 
this  blslgnal  system,  not  only  can  the  pressure  and  temperature  level  be 
detected,  but  also  the  location  of  one  relative  to  the  other  could  be  deter- 
mined. In  addition  to  the  experiments,  analyses  have  been  performed 
during  the  project  to  estimate  the  sources  of  the  measured  temperatures. 
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Considerations  have  been  given  to  effects  such  as  inlet  shear  heating, 
compressive  heating,  and  heating  due  to  disk  slippage. 

The  project  was  conducted  over  a 2-year  period.  The  first  year 
was  devoted  primarily  to  the  development  of  the  bisignal  transducer  and 
preliminary  thermal  analyses.  The  second  year  has  involved  numerous 
temperature  measurements  for  four  selected  lubricants  under  several  con- 
ditions of  speed,  load,  and  disk  slippage.  In  addition  to  the  temperature 

measurements,  traction  measurements  under  slip  conditions  have  been  made  . 

to  aid  in  analyzing  thermal  effects  on  contact  temperatures.  ] 
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EXPERIMENTAL  EVALUATIONS 


The  objective  of  the  experiments  has  been  to  develop  fundamental 
data  to  be  used  for  determining  the  role  of  temperature  in  the  EHD  process. 
IVo  general  types  of  experiments  have  been  conducted  as  follows: 

(1)  Traction-slip  experiments  on  the  type  described 
in  Reference  (1)^^^ 

(2)  Temperature-pressure  measurements  of  the  type 
described  in  References  (2)  and  (3) 

Four  lubricants  were  selected  for  the  experiments: 

(1)  MIL-L-23699 

(2)  XRM-109  (Synthetic  Paraffinic) 

(3)  Commercial  traction  fluid 

(4)  Polyphenyl  ether  (5P4E)  . 

Approximate  properties  for  these  lubricants  are  given  in  Table  1. 

Traction-Slip  Experiments 

The  traction  experiments  were  conducted  using  the  apparatus  shown 
schematically  in  Figure  1.  The  apparatus  consists  of  two  independently- 
driven  disks  which  can  be  loaded  together  into  lubricated  contact.  The 
disks  are  on  stub-shafts  which  fit  into  drive  shafts  mounted  in  precision 
(ABEC-7)  angular  contact  bearings.  The  rotor  for  the  drive  motors  are  inte- 
gral with  the  disk  shafts.  These  motors  are  high-frequency  induction  motors, 
the  speed  of  which  can  be  varied  by  changing  the  input  frequency. 

The  disks  are  jet-lubricated  by  the  test  fluid  at  a preset  tempera- 
ture. Disk  temperature  is  monitored  by  a thermocouple  located  outboard  of 
the  upper  disks.  The  geometry  and  load-pressure  characteristics  of  the  disk 
are  summarized  in  Table  2. 


(a)  References  are  listed  on  Page  41. 
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FIGURE  1.  ROLLING  DISK  APPARATUS,  SET-UP  FOR  X-RAY  FILM  THICKNESS  AND 
TRACTION  MEASUREMENTS 
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TABLE  2.  TYPICAL  PARAMETERS  FOR  TEMPERATURE- 


Parameter 

Symbol 

Metric  Value 

Operating  Condition 

• Load 

W 

450,  1400  N 

• Contact  Pressure 

Ph 

. 7 GPa,  1.1  GPa 

• Speed 

2900  rpm 

Contract  Dimensions 
at  . 7 GPa 

a,b 

233,  1259  pm 

Contract  Dimensions 
at  1.1  GPa 

a ,b 

351,  1890  pm 

Disk  Geometry 

• Radius 

R 

.038  pm 

• Upper  Disk  Crown 
Radius 

R 

.28  m 

TRACTION  STUDIES 

English  Value 

100,  315  lbs 
10^,  1.5  X 10^  psi 
2900  rpm 

.0091,  .0491  inch 
.0137,  .0737  inch 

1.5  inch 

11  inch 
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In  the  traction  slip  experiments,  disk  slippage  is  achieved  by 
reducing  the  input  frequency  to  upper  disk  drive  motor  relative  to  the 
lower  drive  motor.  The  attendant  traction  is  detected  by  a load  cell  con- 
straining the  upper-disk  support  unit.  Slip  is  monitored  using  a differ- 
ential electronic  speed  indicator.  The  output  from  this  indicator  is  fed 
into  the  x-axis  of  an  x-y  recorder  and  the  output  from  the  traction-load 
cell  is  fed  into  the  y-axis.  With  this  arrangement,  a traction-slip  curve 
can  be  generated  for  each  condition  of  interest. 

Traction  slip  data  as  extraced  from  x-y  recordings  are  given  in 
Figures  2 through  8 for  the  four  lubricants  at  several  load  and  temperature 
conditions  using  76  mm  diameter  disks.  As  would  be  expected,  an  increase 
in  load  yields  an  increase  in  the  level  of  the  measured  traction  for  a 
given  slip  condition. 

There  are  two  major  regions  of  interest  in  the  traction  slip 
curves  as  follows  : 

(1)  Near  zero  slip,  the  traction  slip  curves  appear 
to  be  linear,  suggesting  Newtonian  lubricant 
behavior.  For  a Newtonian  lubricant,  the  shear 
stress  (traction)  is  proportional  to  shear  rate 
(slip) . 

(2)  At  higher  values  of  slip,  the  traction  no  longer 
rises  linearly  with  slip.  This  nonlinear  be- 
havior could  be  due  to  either  non-Newtonian  ef- 
fects or  to  thermal  heating. 

It  is  interesting  to  note  that  the  traction  behavior  of  the  XRM-109 
and  the  MIL-L-23699  lubricant  are  somewhat  similar  despite  large  differences 
in  their  base  viscosities.  Also,  the  polyphenyl  ether  is  similar  to  the 
traction  fluid.  However,  the  traction  fluid  exhibits  much  greater  traction 
than  the  MIL-L-23699  even  though  their  viscosities  are  nearly  the  same. 
Traction  is,  of  course,  a complicated  phenomenon  involving  EHD  film  thick- 
ness, contact  zone  rheology,  and  temperature.  One  goal  of  the  project  has 
been  to  evaluate  the  contact  zone  temperature  under  slip  as  well  as  pure 
rolling  conditions.  The  temperature  under  slip  conditions  will  have  a 
profound  impact  on  the  shape  of  the  traction  curve. 
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Temperature-Pressure  Experiments 


Production  of  Thin- 
Film  Transducer 


The  temperature-pressure  (bisignal)  transducer  has  been  described 

(2  3) 

previously  ' and  is  shown  pictorially  in  Figure  9.  This  thin-film  trans- 
ducer is  coated  directly  onto  the  lower  disk  of  the  twin  disk  apparatus. 

This  process  consists  of  the  following  steps: 

(1)  First,  a layer  of  alumina  (^20^)  is  sputtered 
onto  the  disk  surface. 

(2)  A manganin  pressure  transducer  is  vapor-deposited 
onto  the  alumina  layer.  The  shape  of  this  manganin 
coating  (see  Figure  9)  has  been  designed  to  allow 
the  active  region  of  the  transducer  to  be  located 
in  the  axial  center  of  the  disk  surface. 

(3)  A layer  of  alumina  is  sputtered  on  top  of  the 
manganin. 

(4)  The  titanium  thermistor  is  vapor-deposited  onto  the 
last  alumina  layer. 

All  of  the  deposition  processes  are  done  in  a controlled  environ- 
ment broken  only  between  Steps  2 and  3.  The  shape  of  the  manganin  and 
titanium  coatings  is  controlled  with  the  aid  of  masks  manipulated  from  out- 
side the  vacuum  chamber.  The  final  dimensions  of  the  active  portions  of 
the  transducers  are  about  50  x 200  ym  by  .1  ym  thick. 

Calibration  of  Blslgnal 
Transducer 


Theoretically,  as  the  blsignal  transducer  passes  through  the  con- 
junction region  between  the  disks,  the  resistance  change  of  the  titanium  due 
to  temperature  and  of  the  manganin  due  to  pressure  is  detected  simultaneously. 

The  circuitry  used  is  a Wheatstone  bridge  with  a high-speed  oscilloscope 

(3) 

as  the  galvanometer.  However,  as  discussed  the  thermistor  has  been  found 
to  be  somewhat  pressure  sensitive.  Therefore,  a special  pressure  compensating 
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FIGURE  9.  PICTORIAL  DRAWING  OF  P-T  TRANSDUCER 
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circuit  has  been  developed  as  shown  in  Figure  10.  Here,  the  signal  from 
the  pressure  transducer  is  used  to  compensate  for  the  titanium  pressure 
effects.  By  adjusting  the  variable  resistors,  this  pressure  effect  can  be 
min Imlzed. 

The  temperature  transducer  was  calibrated  for  temperature  effects 
by  spraying  oil  at  a known  temperature  onto  the  thermistor.  The  temperature 
calibration  curve  for  the  thermistor  is  given  in  Figure  11;  1 mv  on  the 

oscilloscope  represents  an  8 C temperature  rise.  This  calibration  is  con- 
sistent with  previously  reported  data.  For  the  pressure  calibration,  the 
output  was  compared  with  the  known  loading.  The  calibration  is  approximately 
.18  GPa/mv. 

Bislgnal  Transducer  Output 

A typical  oscilloscope  trace,  obtained  with  the  bisignal  transducer, 
is  given  in  Figure  12  for  a MIL-L-23699  lubricant  at  a loading  of  1.1  GPa 
contact  pressure.  Since  the  compensating  c:'rcuit  (Figure  10)  precludes 
recording  pressure  and  temperature  at  exactly  the  same  time,  a double  ex- 
posure photograph  has  been  made;  that  is,  an  image  of  the  temperature  trace 
immediately  followed  by  the  pressure  trace  is  taken  on  a single  photograph. 

A single  magnetic  pickup  is  used  to  trigger  the  oscilloscope  and,  hence, 
to  provide  horizontal  location  of  the  pressure  or  temperature  trace  on  the 
scope  screen.  The  shutter  speed  of  the  camera  was  slow  enough  that  several 
traces  of  the  pressure  and  temperature  distributions  are  shown  in  the  figures. 

The  inlet  region  in  all  of  the  photographs  is  on  the  left  side. 

The  temperature  grows  in  the  inlet  region  and,  for  pure  rolling,  peaks  near 
the  peak  pressure  gradient.  The  shape  of  the  pressure  trace  is  very  near  a 
semlelllptic  shape  which  is  similar  to  the  Hertzian  theory  for  dry  contact. 

Near  the  exit  region,  a "dip"  occurs  in  the  temperature  curve. 

This  is  thought  to  be  a result  of  the  pressure  correction  circuit  and  a 
slight  misalignment  of  the  temperature  transducer  relative  to  the  pressure. 

A misalignment  of  30-50  pm  could  cause  a small  region  where  the  pressure 
compensation  would  cease  to  exist  and  cause  an  artificial  drop  in  the 
temperature  trace.  Outside  the  pressure  region,  the  thermistor  returns 
to  the  proper  level  since  pressure  effects  are  negligible. 


Variable 

resistor 


FIGURE  10.  SELF  PRESSURE-COMPENSATING  CIRCUIT  FOR 
TEMPERATURE  EXPERIMENTS 


mV  Output  on  Scope 
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FIGURE  11.  TEMPERATURE  CALIBRATION  CURVE 
FOR  THERMISTOR 
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Pressure  Curve 


Temperature  Curve 


Region  of  Uncorrected 
Pressure  Effects  on 
Temperature  Transducer 


FIGURE  12. 


TYPICAL  TEMPERATURE- PRES SURE  SCOPE  TRACE 
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Slip  Effects  on 
Measured  Temperature 


Temperature-pressure  traces  are  presented  in  Figures  13  through 
16  for  the  four  selected  lubricants  (see  Table  1)  for  one  load  (400  N) , 
one  speed  (2800  rpm) , and  ambient  temperature  at  various  slip  conditions 
including  pure  rolling.  Figure  13  presents  the  temperature  data  for  a MIL- 
L-23699 lubricant.  Under  pure  rolling  conditions,  the  temperature  rises 
to  a peak  of  32  C near  the  point  of  maximum  inlet  pressure  gradient  and 
is  followed  by  a gradual  temperature  drop.  The  growth  of  temperature  in 
the  inlet  region  is  caused  by  a combination  of  rolling  shear  on  the  oil 
and  compressive  heating.  Figures  13(b)  to  (d)  show  that  slip  does  not  af- 
fect the  inlet  zone  temperature  but  does  materially  effect  the  temperature 
at  the  exit  region.  For  example,  at  17  percent  slip,  the  exit  temperature 
has  increased  by  30  C.  The  temperature  data  for  the  XE<M-109  (Figure  14) 
are  similar  to  the  MIL-L-23699.  For  the  traction  fluid  (Figure  15)  and 
polyphenyl  ether  (Figure  16) , the  temperature  rise  under  pure  rolling  is 
about  the  same  as  the  other  two  lubricants.  However,  the  temperatures  do 
not  decay  through  the  contact  region  for  these  fluids  even  under  pure 
rolling  conditions.  An  enormous  temperature  rise  occurs  for  the  traction 
fluid  under  high  slip  [Figure  15(a)  ] due  to  the  high  level  of  traction 
associated  with  this  fluid  (see  Figure  5) . 

Figure  17  is  a plot  of  exit  temperature  rise  as  a function  of 
slip  for  the  four  lubricants.  As  can  be  observed,  the  MIL-L-23699  and 
the  XRM-107  are  quite  similar.  Also,  it  can  be  seen  that  the  temperature 
rise  with  the  traction  fluid  is  significantly  larger  than  for  the  other 
lubricants . 

In  general,  the  high  magnitude  of  the  temperature  rise  due  to 
slip  could  explain  much  of  the  nonlinearities  seen  in  traction-slip  curves 
at  high  slip  conditions. 

Rolling  Speed  Effects 

Figure  18  presents  temperature-pressure  traces  for  a rolling  speed 
of  1400  rpm  (l.e.  , approximately  half  the  speed  of  the  previous  data).  This 
reduction  in  speed  significantly  reduced  the  temperature  rise  under  pure 
rolling  conditions.  Figure  19  is  a bar  graph  showing  speed  effects  on 
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(c)  7%  Slip  (d)  17%  Slip 


FIGURE  13.  TEMPERATURE-PRESSURE  DATA  FOR  A MIL-L-23699  LUBRICANT 

Speed  = 2800  rpm  (.076  m disk) 

Temperature  Calibration  = 16  C/div 
Load  = . 7 GPa 
Oil  Temperature  = 35  C 
Scope  Sweep  = 20  ps/div 
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(a)  "Pure"  Rolling  Disk 


(b)  7%  Slip 


(c)  15%  Slip  (d)  23%  Slip 


FIGURE  14.  TEMPERATURE  PRESSURE  DATA  FOR  XRM-109  LUBRICANT 

Speed  = 2800  rpm  (.076  m disk) 

Temperature  Calibration  = 16  C/div 
Load  = . 7 GPa 
Oil  Temperature  = 45  C 
Scope  Sweep  = 20  ps/div 
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(c)  2.2%  Slip  (d)  13%  Slip 


FIGURE  15.  TEMPERATURE  PRESSURE  DATA  FOR  TRACTION  FLUID 

Speed  = 2000  rpm  (.076  m disk) 

Temperature  Calibration  = 16  C/div 
Load  = . 7 GPa 
Oil  Temperature  = 35  C 
Scope  Sweep  = 20  ps/div 
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(c)  lOX  Slip 


(d)  13%  Slip 


FIGURF.  16.  TFMPI-.RATURE  PRESSURE  DATA  FOR  POl.YPHF.NYL  ETHER  (5P4E) 

Spri’d  = 2800  rpm  (.076  m disk) 

Ti-nipcraLure  Calibration  = 16  C/div 
Load  . 7 CPa 
Oil  ToniporaLure  s 45  C 
Scope  Sweep  = 20  ps/div 


(a)  Traction  Fluid 


(b)  MIL-L-23699 


FIGURE  18.  TEMPERATURE  PRESSURE  DATA 
FOR  SELECTED  LUBRICANTS 

Speed  = 1400  rpm  (.076  m disk) 
Temperature  Calibration  = 16  C/div 
Load  = . 7 GPa 
Scope  Sweep  = 40  ps/div 
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FIGURE  19.  EFFECT  OF  SPEED  ON  "PURE"  ROLLING  TEMPERATURE  RISE 
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Inlet  temperature  rise.  The  temperature  drops  roughly  as  the  square  root 
of  the  speed. 

Load  Effects 


Figure  20  presents  temperature-pressure  traces  for  a load  of  1.1 

GPa.  For  this  load,  the  temperature  appears  to  rise  to  about  50  C under 

pure  rolling  for  most  of  the  lubricants,  as  opposed  to  30  C for  the  .7  GPa 

load.  It  is  possible  that  this  load  effect  on  temperature  could,  at  least, 

partially  explain  anomolous  load  effect  on  EHD  film  thickness  seen  in  various 
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experiments  such  as  BCL-X-ray  experiments. 


(a)  MlL-L-23699 


(b)  XRM-109 


(c)  Traction  K]iiid 


FIGURE  20.  TEMI’ERAniRE-l’RESSURE  DATA  FOR  FOUR  SKl.F.CTED  LUBRICANTS 

Speed  = 2H00  rpm  (.076  m disks) 

Temperature  Calibration  = 16  C/div 

Load  =1.1  Cl'a 

St'opt^  Swi'cp  = 20  ps/div 
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THEORETICAL  EVALUATIONS 


The  purpose  of  the  analyses  has  been  to  understand  the  nature 
of  the  EHD  temperature  and  to  assess  the  validity  of  the  temperature 
measurements.  As  shovm  in  the  experimental  section,  these  measurements 
have  indicated  that  very  high  temperature  rises  can  occur  under  pure  roll- 
ing as  well  as  rolling/sliding  EHD  contacts.  Temperature  rises  on  the  order 
of  30-150  C can  have  a profound  effect  as  lubricant  rheology  and,  hence,  on 
tractions  and  EHD  film  thickness. 

Rheology  Evaluations  for  Tractions 

Very  recently,  rheology  has  been  the  topic  of  various  researchers, 
including  Johnson^^\  Winer^^\  Miller^^\  and  Montrose^^\  As  a result 
of  the  rese_.ch  by  these  scientists,  very  interesting  non-Newtonian,  limit- 
ing shear,  transition,  time  transient,  and  viscoelastic  effects  have  been 
Isolated.  The  purpose  of  the  BCL  study  has  been  to  establish  baseline 
data  to  be  used  to  assess  the  effect  of  temperature  on  lubricant  behavior. 
More  work  is  needed  in  the  area  of  interrelating  these  rheology  analyses 
with  thermal  effects. 

Under  low  slip  condtions,  shear  heating  due  to  slip  should  not 
be  a major  factor  in  controlling  EHD  tractions.  By  analyzing  the  slope 
of  the  traction-slip  curves  at  zero  slip,  then,  it  should  be  possible  to 
determine  a pseudo-viscosity  for  the  lubricant  under  rolling  conditions. 

This  viscosity  is,  of  course,  related  to  the  inlet  temperature  rise  and 
could  be  related  to,  say,  glass-transition  effects.  For  these  evaluations, 
we  have  chosen  to  simply  label  this  as  a limiting  viscosity  ^or  specific 
operating  conditions. 

To  calculate  this  limiting  viscosity,  it  .is  helpful  to  assume 
that  the  tractive  force,  F^,  between  the  disk  can  be  expressed 

"t-  • <« 

Here  t Is  shear  stress,  A is  contact  area,  AV  is  slip,  h is  the  EHD  film 
thickness,  and  p(P)  is  the  viscosity  of  the  lubricant  at  the  contact  pressure 
conditions.  If  p is  assumed  to  be  a limiting  viscosity,  then 
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W 


AAV 


(2) 


wliore,  for  elliptical  contacts, 


A = irab  . (3) 

Here,  b and  a are  the  major  and  minor  axes  of  the  Hertz  contact  ellipse. 

Values  of  p,  as  computed  using  Equation  (2)  are  given  in  Table  1.  The  film 
tliickness  was  computed  using  a semiemperical  equation  developed  at 

Temperature  Rise 
Due  to  Sliding 

In  the  first  summary  report,  a theoretical  analysis  of  the  tem- 
perature rise  due  to  slip  was  developed  (see  Appendix  A).  This  analysis 
was  based  on  heat  generation  due  to  slip  (only)  and  heat  loss  due  to  conduc- 
tion across  the  lubricant  film.  The  temperature  rise  (T)  could  be  written 
in  the  form 


T = 3.1 


/S 


y (p)a 


K 6Vpc 
s 


h 


(4) 


where  K and  K are  the  lubricant  and  steel  conductivity,  6 is  the  tempera- 

Lt  S 

ture  coefficient  of  viscosity,  p is  the  density  of  the  disk,  c is  the  speci- 
fic heat  of  the  disk,  V is  the  surface  velocity,  a is  the  half  width  of 
contact,  and  h is  film  thickness.  Reasonable  values  for  these  parameters 
are 


= .30  N/sec  C 

K = 35  N/sec  C 
s 

6 = .03/C 

V = 11  m/sec 
pc  = 3.5  X 10^  N/m^C 
= 230  pm. 


a 
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Then, 


T = 4,0  X 10“6  Ax  ^ C 

n N 


(5) 


Equation  (5)  can  be  used  to  compute  the  temperature  rise  due  to 
sliding  contact.  The  results  of  these  estimates  are  summarized  in  Table  3 
in  comparison  with  the  data  taken  from  Figure  17.  As  can  be  observed,  good 
agreement  exists  between  theory  and  experiment. 


Temperature  Rise 
Due  to  Rolling 


The  energy  equation  of  the  lubricant  film  can  be  written 


(9) 


„ 8T  ^ 32t 

9x  ~ ^ 2 

3y 


Q + Q 
s c 


(6) 


where  and  c^^  are  the  density  and  specific  heat  of  the  lubricant, 
is  the  thermal  conductivity,  x is  the  tangential  coordinate,  y is  normal 
coordinate,  V is  velocity,  is  the  heating  due  to  shear,  and  is  the 
heating  due  to  compression. 

It  is  helpful  to  scale  the  left  side  of  the  equation  in  the 
following  way; 


^ _ ^^o  B^T  o + n 

a 3C  h ^ Tn^  " ^s  ^c  • 
o 

Typical  values  for  the  constants  are  as  follows: 

pc  = 1.5  X 10®  N/m^C 

a = 230  pm 
V = 11  m/sec 
= .3  N/sec  C 

h = 5 pm  (max)  . 


(7) 
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TABLE  3.  ESTIMATE  OF  SURFACE  TEMPERATURE  RISE  DUE  TO  SLIDING 


Lubricant 

Viscosity 

at  Pressure 

Film  Thickness, 
10“'’  m 

Temperature  Rise  C 
at  20%  Slip 

N Sec 

Cp 

Calculated 

Measured 

MIL-L-23699 

lb 

.16  X 10^ 

1.35 

27 

30 

XRM  109 

100 

I.O  X 10^ 

5.1 

18 

25 

Traction  Fluid 

320 

3.2  X 10^ 

1.46 

108 

130 

Polyphenyl  Ether 

2330 

23  X 10^ 

4.  7 

90 

75 
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Using  these  values  then, 


7 X 10 


10 


r ^ 

o 6C 


- 1.2 


10 

X 10  T 


o 9n 


3^-T  _ „ 
7. — Q, 


+ Q 
s 


(8) 


In  the  inlet  zone  (where  h > h^)  for  at  least  two  of  the  lubricants 
evaluated,  an  approximation  to  Equation  (5)  can  be  written 


pcV 


3T 

9x 


Q + Q 


(9) 


Now,  it  can  be  shown  that 


T “ "1^  y (pure  rolling  equilibrium) 

dp  *^“^0 

= 12  p(p)  u (Reynolds  equation) 

t2 

Q = — (shear  heating) 
s p ^ 

(Burton* s^^^^  equation) 

h = + «,n  (5  + 

(Hertz  inlet  shape) 

Here, 


(10) 

(11) 

(12) 

(13) 

(14) 


p = pressure 

h = film  thickness  at  any  point  in  the  inlet  region 

T = shear  stress  in  the  lubricant  film 

B = a compressibility  constant 

T = absolute  temperature  (BT  « .1) 

3 A 


a = half-width  of  constant 
C = x/a 

R = the  relative  radius  of  the  disk 


p = viscosity  of  the  lubricant  [p  = p^  exp(Yp-6T)l  . 
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form: 


Equations  (9)  through  (14)  can  be  combined  to  the  following 


and 


where 


. , 9n  V6  (h-h  )2  12  p VBT,6  h-h 

^ _ __o o o A o 


li  - 12tJVy  *^~^o 
9C  ~ ~ Ip  h'^ 


^ -YP  , , 6T 
(P  = e and  = e 


(15) 


(16) 


(17) 


Equations  (14)  and  15)  are  amenable  to  solution  by  numerical  integration 
such  as  a Runga  Kucta  scheme. 

Temperature  predictions  using  Equation  (14)  through  (16)  are 
shown  in  Figure  21.  Except  for  the  MlL-L-23699  evaluation,  an  insta- 
bility in  the  temperature  solutions  occurred  as  a result  of  very  large 
pressure  build-up.  The  unfortunate  aspect  of  the  equations  are  that 
they  represent  only  the  inlet  region  and  do  not  merge  easily  with  contact 
zone  equations.  The  difficulty  is  that  the  pressure  tends  to  infinite 
which  in  reality  would  deflect  the  surface  and  abate  the  pressure.  The 
magnitudes  of  the  temperature  predictions  in  Figure  (21)  are  higher  than 
measured  since  the  equations  represent  essentially  an  upper  bound  for 
temperature.  However,  the  magnitudes  of  the  predicted  temperatures  illus- 
trate that  the  experimental  data  are  reasonable  and  that  very  high  tempera- 
tures occur  in  EHD  contacts. 

In  the  contact  region,  compressive  heating  is  the  major  heat 
source  under  pure  rolling  conditions.  If  we  Ignore  the  shear  stress 
effect  in  Equation  (5),  then  it  can  be  seen  that 


If,  for  example,  BT^  = .1  and  p = .7  x 10^  N/m^  (.7  GPa) , then,  T = 46C 
which  would  be  the  adiabatic  temperature  rise  due  to  compression  in  the 
contact  region.  This  estimate  further  illustrates  that  the  measured 
temperatures  were  reasonable. 
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FIGURE  21.  PREDICTED  INLET  ZONE  TEMPERATURE 
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Critique  of  Techniques 

The  observations  made  with  the  temperature  transducer  are  very 
significant  with  regard  to  lubricant  behavior  in  concentrated  contacts  and 
on  the  whole  EHD  process.  It  is  important,  then,  to  challenge  the  technique 
itself  to  assess  its  reliability.  The  temperature  calibration  of  the 
transducer  is  relatively  straightforward  and  seems  to  be  consistent  for 
a range  of  temperatures.  Further,  such  limitations  as  time  response  of 
the  transducers  can  easily  be  discarded  since  the  pressure  transducer  (with 
a similar  circuit)  is  responding  to  faster  events  than  the  thermistor  with 
no  obvious  limitations. 

Pressure  effects  on  the  titanium  have  been  examined  and  efforts 
to  compensate  for  these  effects  have  been  made.  In  general,  for  the  .7 
CPa  experiments,  we  are  confident  that  good  compensation  has  been  made. 

At  higher  loads,  it  appears  that  the  pressure  compensation  is  slightly 
insufficient  and  the  temperatures  are  actually  higher  than  measured.  The 
pressure  effect  on  the  thermistor  was  checked  by  inserting  thin  capacitance 
paper  between  the  disks  and  rolling  them  together  under  a 400  N load.  With 
the  compensating  circuit,  the  titanium  pressure  effect  was  only  about  +15 
percent  of  the  manganin  pressure  effect.  In  all  the  experiments  reported 
here,  the  thermistor  output  was  taken  on  the  2 mv/dlv  scale  of  the  oscillo- 
scope and  the  pressure  output  was  taken  on  the  1 mv/div  scale.  If  the 
level  of  the  pressure  and  temperature  signal  were  the  same,  the  pressure 
error  on  the  temperature  measurement  would  be  only  about  7 percent  of  the 
temperature  signal. 

A final  point  is  that  the  theoretical  temperature  rise  due  to 
slip  (for  steel  disks)  is  consistent  with  the  measurements.  This  further 
tends  to  corroborate  both  the  calibration  of  the  transducer  and  the  simili- 
tude between  the  disk  (with  the  alumina  layer  on  it)  and  pure  steel  disks. 
Since  the  published  thermal  conductivity  for  alumina  is  very  similar  to 
steel,  this  is  quite  reasonable. 
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CONCLUSIONS 


The  objective  of  the  project  has  been  to  evaluate  the  role  of 
temperature  in  the  EHD  process.  Based  on  the  project  experiments,  the 
following  conclusions  were  made. 

(1)  The  inlet  temperature  rise  was  approximately  30  C 
for  all  four  lubricants  evaluated  at  .7  GPa. 

(2)  The  temperature  drop  in  the  contact  region  (under 
pure  rolling  conditions)  depends  heavily  on  the 
lubricant  and  apparently  the  lubricant  film  thick- 
ness (l.e.,  the  volume  of  heated  fluid). 

(3)  The  temperature  rise  under  pure  rolling  varies 
roughly  as  the  square  root  of  rolling  speed  and 
increases  with  increasing  load. 

(4)  Contact  zone  temperatures  as  high  as  200  C have 
been  measured  under  high  slip  conditions  with  a 
traction  fluid. 

Analyses  in  connection  with  the  experiments  indicate  that 

(1)  The  pure-rolling  temperature  rise  can  be  explained 
based  on  inlet  shear  and  compression  heating. 

(2)  The  measured  temperature  rise  due  to  slip  is  reason- 
able and  consistent.  Further,  this  temperature 
rise  could  be  used  to  explain  much  (and,  perhaps 
all)  of  the  drop  off  in  traction  due  to  slip, 
especially  for  thick  film  lubrication. 

(3)  The  observed  drop  off  in  temperature  with  decreasing 
speed  makes  the  results  reasonable  in  comparison 
with  other  researchers. 

(4)  The  increase  in  temperature  with  increasing  load 
hints  that  temperature  is  a major  reason  for  anamolous 
decreasing  film  thickness  with  load,  as  observed  with 
X-ray  measurements. 
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(5)  The  overall  magnitude  of  the  temperature  measure- 
ment makes  isothermal  EHD  film  thickness  theories 
very  suspect  for  real  bearing  contacts 

The  results  of  this  project  point  to  several  areas  where  more 
research  would  be  fruitful.  These  areas  include 

(1)  An  analysis  of  the  temperature  rise  in  pure  rolling 
for  a range  of  loads  and  speeds  to  evaluate 

(a)  Effect  of  load/temperature  on  film  thickness 

(b)  Effect  of  speed/temperature  of  film  thickness, 

(2)  Experiments  to  determine  the  effect  of  non-Newtonian 
lubricants  such  as  polymer  thickened  oils  or  greases 
on  contact  zone  temperatures  and  bearing-gear 
performance , 

(3)  A better  definition  (both  analytical  and  experimental) 
of  the  role  of  temperature  in  traction  to  aid  in  de- 
signing traction  drive  units. 
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Energy  Equation  for  the  Lubricant  Film 

The  purpose  of  this  analysis  is  to  develop  a simple  method  for 
estimating  surface  temperature  rise  due  to  EHD  tractions.  The  approach 
presented  here  is  similar  to  the  one  developed  by  Kannel  and  Bell.^^^^ 
However,  in  that  paper,  a lubricant  model  using  an  exponential  temperature 
relationship  was  used.  This  model  resulted  in  an  overly  complicated  expres- 
sion for  surface  temperature  distribution.  A more  convenient  lubriant  model 
is  of  the  form 


where 


p^exp(YP) 

^ " 1+8 (T-T^)  ’ 

U “ the  lubricant  viscosity 
Vi^=  viscosity  in  the  inlet  region 
Y = pressure  viscosity  exponent 
p = pressure 

6 = temperature  coefficient  of  viscosity 

T = temperature 

Tj*  inlet  zone  temperature  . 


(A-1) 


Ignoring  convection  of  heat  along  the  film  in  favor  of  conduction 
across  the  film,  the  energy  equation  for  a Newtonian  lubricant  model  can 
be  wrltten^^^^ 


^ 2 
a^T  . / 3u\  . 


(A-2) 


A-2 


where , 

T = temperature 

Kj^=  thermal  conductivity  of  lubricant. 


It  is  helpful  (though  not  mandatory)  to  assume  that  the  temperature  distri- 
bution in  the  lubricant  film  is  symmetrical  about  the  centerline.  Also, 
(12) 

it  is  assumed  that  the  shear  stress  is  constant  across  the  film.  With 
these  assumptions  and  the  definition  of  a thermal  loading  parameter,  T , 

Ij 

given  by 


/u . 6 

(U2-U1)  , 


(A-3) 


there  results  from  Equation  (A-2) 


or 


Here , 


By 


s 


arccot 


2(1+6T[0])~ 

Tj^exp(Yp/2)  ’ 


By 


s 


2h6 


exp 


IE 

2 • 


(A-4) 


(A- 5) 


•Uj^  = sliding  velocity 

h = film  thickness 

= thermal  conductivity  of  lubricant 

s = surface  temperature  or  temperature  gradient. 


Heat  Transfer  to  the  Solid 


Using  the  equation  of  Cheng 
solid  can  be  written 


(13) 


the  heat  transfer  to  the 


T -T 
s 


i 


/ U dC 

-a  s /x-^ 


(A-6) 


Combining  (A-5)  and  (A-6) , the  following  expression  for  surface  temperature 


results : 


where  , 


r . J-  yS  ? expriEia)  _ss_ 


k = diffuslvity  of  steel 


“ thermal  conductivity  of  steel 
p * density  of  the  lubricant 
c • specific  heat 

a = half-width  of  the  contact  region 


= inlet  temperature  . 


Approximate  Solution  for 
Surface  Temperature 


The  pressure  distribution  between  heavily  loaded  EHD  contacts 
can  be  approximated  to  a high  degree  of  accuracy  by  the  equation 


(A- 7) 


p-Ph/i-  i • 


(A- 8) 


If  we  assume  that  heating  occurs  very  near  the  center  of  contact,  then. 


1--  Mexp 


^Ph  r 1 ^Ph  /x\n 

— [-  2 ~ [l)  J- 


Now,  by  defining 


»-(^)  (f-f) 


(A-9) 


(A-10) 


there  results 


T -T  a — — 

s i <Sh  K 

s 


where , 


^ exp  ^(8-0^)  ^d0  . 


f(6)  , 


(A-11) 


(A-12) 


Valui's  of  f(3)  are  presented  in  Figure  A-1.  The  peak  value  is  1.3. 
Combining  Equation  A-3  and  A-11,  there  results 


T =21  e^^ha  ^ / 2 V 

fiVpc  h Vyp.  / 


,1/4 


max 


(A-13) 


For  an  elliptical  contact  region,  it  can  be  shown  that 


"ave  = X i • 


or* 


(A-14) 


Thus,  (approximately) 


P = — / p e^P  (-2Tiab  pdp)  . 
ave  ab  o o i i / 


Ph  ^"Ph^' 

Pq  ^ ' Pave  ^(yp,-!)  ’ 


(A-15) 


(A-16) 


if  we  assume  that  yPj^  = 3 for  high  pressure  condition.  Then, 


p e^Ph  = 2.25  p 
o 


ave 


(A-17) 


and  Equation  (A-13)  becomes 


T =3 
s 

max 


. A.  Pave' 

6Vpc 


6Vpc  h 


(A-18) 


* See  Reference  (1),  Equation  (21). 


